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Abstract. A density-dependent two-nucleon potential has been derived in the formalism of 
correlated basis function. The effects of 3-particle interactions has been included by integrating 
out the degrees of freedom of the third nucleon. The potential can be easily employed in nuclear 
matter calculations. It yields results in agreement with those obtained from the underlying 
three-body potential. The use of the density dependent potential allowed us to study the effects 
of three-nucleon interactions in symmetric nuclear matter within the Auxiliary Field Diffusion 
Monte Carlo (AFDMC) computational scheme. 



1. Introduction 

The Auxiliary Field Diffusion Monte Carlo (AFDMC) [T] results of Ref. [5] do not lead to a 
decrease of the binding energy predicted by Fermi-Hyper-Netted-Chain (FHNC) and Brueckner- 
Hartree-Fock (BHF) calculations [4J. This suggested that the overstimation of the SNM binding 
energy resulting from FHNC calculation with UIX three body potential [5], may not be solved 
with a better choice of the variational wave function. 

More sofisticated three body potentials, while providing a quantitative account of the energies 
of the ground and low-lying excited states of nuclei with ^4 < 12 [6], give very different results 
in pure neutron matter (PNM) [3]. 

Lagaris and Pandharipande [7] and Friedman and Pandharipande [8] constructed an effective, 
density-dependent, two-nucleon potential (TNI) to take into account three- and many-nucleon 
forces. However, they adopted a purely phenomenological procedure, as the parameters of the 
density dependent functions appearing in this potential were determined through a fit of the 
saturation density, the binding energy per nucleon and the compressibility of symmetric nuclear 
matter (SNM). In general, fitting the potential to results obtained with many-body techniques 
makes the potential itself affected by the approximations of the many body calculations. 



Our potential [9] improves the TNI model since it has been derived from a realistic microscopic 
three-nucleon force, the UIX potential, which provides a good description of the properties of 
light nuclei, although very recently it has been shown [10] that UIX does not reproduce both the 
nd scattering lenght and the vector polarization observables Ay and iTn. However, our method 



can be applied to more refined three nucleon interactions, like the chiral NNLOL described in 

m- 

To obtain the density dependent potential, the average on the degrees of freedom of the 
third particle has been carried out using a formalism suitable to account for the full complexity 
of nuclear dynamics. Our results show that, in doing such reduction, of great importance is 
the proper inclusion of both dynamical and statistical NN correlations, whose effects on many 
nuclear observables have been found to be large jlll I12| . 

We have used CBF and the Fantoni-Rosati (FR) cluster expansion formalism[T3] to perform 
the calculation of the linear- in-density terms of the effective potential, arising from the irreducible 
three-nucleon interactions modeled by the UIX potential. 

The effective potential has been implemented in the AFDMC computational scheme to obtain 
the EoS of SNM, while similar calculations using the UIX potential are not yet possible, due 
to the complexities arising from the commutator term. In addition, the density-dependent 
potential can be used to include the effects of three-nucleon interactions in the calculation of the 
nucleon-nucleon scattering cross section in the nuclear medium. The knowledge of this quantity 
is required to obtain a number of nuclear matter properties of astrophysical interest, ranging 
from the transport coefficients to the neutrino emission rates [141 ES] • 

2. Two- and three- nucleon forces 

The Argonne vg [16j two-body potential model is given by 

% = E^^(-..)of, (1) 
p=i 

where 

Of~^"^ = {l,(7ij,Sij,Lij ■ Sij) (g) (l,Tjj) . (2) 

In the above equation, dij = cTi ■ cTj and Tij = Ti ■ tj, where (Xj and Tj are Pauli matrices acting 
on the spin or isospin of the i-th particle, while Sij is the tensor operator, Ljj is the relative 
angular momentum and Sij is the total spin of the pair. 

We used the so called Argonne Vg and Argonne v'q potentials, which are not simple truncations 
of the Argonne uig potential [TT], but rather reprojections [T8] obtained by refitting the scattering 
data and the deuteron binding energies. In all light nuclei and nuclear matter calculations the 
results obtained with the v'g are very close to those obtained with the full vis- 

Using a nuclear Hamiltonian including only two-nucleon interactions leads to the 
underbinding of light nuclei and overestimating the equilibrium density of nuclear matter. Hence, 
the contribution of three-nucleon interactions must be taken into account, by adding to the 
Hamiltonian the corresponding potential. One of the most widely used is Urbana IX (UIX) [5] 
potential, consisting of two terms: the Fujita and Miyazawa |19j attractive two-pion exchange 
interaction V'^'^ and the purely phenomenological repulsive term 

y27r ^ ^27r ^ (^{Xij , Xjk}{Tij ,Tjk} + -[Xij , Xjk][Tij , Tjk]^ 
cyclic 

V^ = UoY, T\m^r,j)T\m^rjk) ■ (3) 

cyclic 

The spin structure of V'^'^ is given by 

Xij = Y{mT,r)aij + T{m^r)Sij . (4) 



The radial functions read 



Y{x) = —iY{x) , T{x) = (l + - + \)Y{x)iT{x) (5) 

X \ X X'^ / 

where ^y(x) = ^t{x) = 1 — exp(— cx^) are short-range cutoff functions. The parameters ^27r and 
C/q are varied to fit the observed binding energies of and ^He and to reproduce the empirical 
nuclear matter saturation density, while the cutoff parameter is kept fixed at c = 2.1 fm~^. 



3. Correlated Basis Theory and Cluster Expansion technique 

In the correlated basis theories of Fermi liquids [201 121j . the expectation value of the two-body 
potential can be written in the form 

{v) = \p^ j dfidf2vf2 9u ' (6) 



where 



p 



p _AiA- 1) TVi2 / dx3 . . . dxAnF^Ol^F<P, 



pi JdX^F^F^o ' 

are the operatorial components of the two-body distribution function. The uncorrelated 
wavefunction $o in nuclear matter is conveniently chosen to be a Slater determinant of plane 
waves. 

The structure of the correlation operator F reflects the complexity of the Argonne v'q nucleon- 
nucleon potential |22j : 

A 6 
F = 5 n ^ith 4- = ^ r(r,,)Of^. . (8) 

j>i=l p=l 

The radial functions f^irij), appearing in the definition of the correlation operator are 
determined by the minimization of the energy expectation value Ey = (^'ol^^l^o); which 
provides an upper bound to the true ground state energy Eq. As explained in detail in |9|, 
the calculation of Ey in CBF theories is carried out by i) expanding {^o\H\^q) in powers of 
dynamical correlations h{rij) = f^{rij)'^ — 1, 2f^{rij)f'^{rij), {rij) f^^^ {rij) that vanish in 
uncorrelated matter and ii) summing up the main series of the resulting cluster terms by solving 
a set of coupled integral equations. To accomplish the first of these two steps, an extension of 
the FR cluster expansion [T3], that is able to treat scalar correlations only, was introduced in 
[22j to deal with spin-isospin dependent correlation operators, like those of Eq. dS]). 

The cluster terms are most conveniently represented by diagrams (22]. consisting of dots 
(vertices) connected by different kinds of correlation lines. In particular hij is usually represented 
by a dashed line, '^ffjlfj by a single wavy line, and ffj'^ ffj^^ by a doubly wavy line. In addition 
to the dynamical correlation lines, there are also statistical correlation lines, represented by 
solid oriented lines forming closed loops that never touch each other. They arise from the Slater 
determinant and are associated with the Slater function iij 

Open dots of the diagrams represent the active (or interacting) particles (1 and 2), while black 
dots are associated with passive particles, i.e. those in the medium. Integration over the 
coordinates of the passive particles leads to the appearance of a factor p. 
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Figure 1. Diagrammatic representation of Eq. (|13|) . 



It can be shown that all the diagrams contributing to Ey, hence to (v), are linked. Those built 
with scalar passive bonds only, with the only exception of the so called elementary diagrams, 
can be summed up in closed form by solving the FHNC equations [13j. 

On the other hand, diagrams having one or more passive operatorial bonds are calculated at 
leading order only. This implies that at most two operatorial passive bonds can be attached to 
any internal point, thus only diagrams with Single Operator Chain (SOC) are considered. Such 
an approximation is justified by the observation that operatorial correlations are much weaker 
than the scalar ones. 



4. Density dependent potential 

As for the two-body potential vg, it is possible to write the UIX three body potential in the 
following way 

^123 = E^l'30?23- (10) 

P 

The expectation value of V123 then reads 

^-^ = j,p'E[^^^^dn,vr,,gl,,, (11) 
p 

with 

p A\ Tri23 J rfx4 . . . (ix^<l>|)FtOP,3F$o 

'^^^"(^-3)! p^JdX<^>lF^F<l>o ■ ^ ^ 

We require that the expectation values of V123 and of fi2(p) be the same, implying in turn 

drsVus gus = Y] v^ip) 9i2 ■ (13) 

A diagrammatic representation of the above equation, which should be regarded as the definition 
of the vi2{p), is shown in Fig. [TJ The graph on the left-hand side represents the three-body 
potential times the three-body correlation function, integrated over the coordinates of particle 
3. Correlation and exchange lines are schematically depicted with a line having a bubble in 
the middle, while the thick solid lines represent the three-body potential. The diagram in the 
right-hand side represents the density-dependent two-body potential, dressed with the two-body 
distribution function. Obviously, has to include not only the three-body potential, but also 
the effects of correlation and exchange lines. 



Figure 2. Diagrams contributing to the density-dependent potential. The dashed hnes with 
diamonds represent the first order approximation to g^o^^ifij)-, discussed in the text. 



In the construction of the density dependent interaction we have found that the most relevant 
diagrams are those depicted in Fig. [2] that indeed include both statistical and dynamical 
correlations. To simplify the pictures, the three-body potential acting on particles 1, 2 and 
3 is not explicitely depicted. In order to include higher order cluster terms, we have replaced the 
scalar correlation line /j?-^ with the Next to Leading Order (NLO) approximation to the bosonic 
two-body correlation function: 



c 2 



..NLO, 
ybose 



(rij) = /, 



c 2 



0-/ 



df3hi3h23 



(14) 



The inclusion of both statistical and dynamical correlations plays a fondamental role in 
the determination of a realistic density dependent interaction, as can be seen from Fig [3l 
where different density dependent potentials contributions are plotted and compared with the 
full calculation with the genuine three- nucleon interaction. In particular ?Ji2(/o), obtained by 
averaging over the third particle without statistical and dynamical correlations is far from the 
solid curve of UIX. When only statistical correlations are considered, the relative potential, 
ffK/o), comes closer to the UIX curve. Only when the full set of diagrams of Fig. [2] is included 
in the calculation 



9b, 



^"°(n3)5L"e°(r23)(l 



2Pi3^^3)+45L'f (n3)/c(r23)/(r23) , (15) 
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Figure 3. Contributions of the density-dependent potential to the energy per particle of SNM 
(a) and PNM (b), compared to the expectation value of the three-body potential UIX: {Vi23)/A. 



where /(r23) denotes the sum of non central correlations, the contribution of the density 
dependent potential comes very close to that of the original UIX. A proper treatment of the 
statistical correlation is crucial for the construction of a density dependent potential to be treated 
as an additive term to the standard two-body potential. Some care is needed for the density 
dependent potential to reproduce the exchange loop involving particles 1, 2 and 3 with the 
appropriate symmetry factor (see j9] for a thorough discussion). 

5. Numerical results: PNM and SNM equation of state 

As explained in [9], the correlation functions are determined by minimizing the two-body cluster 
term of the energy expectation value for a set of parameters dc, dt, (3p and Op. The energy 
expectation value Ey calculated in full FHNC/SOC approximation has then to be minimized 
with respect to the variation of these parameters. To this aim a Simulated annealing [23] 
optimization algorithm has been implemented [9]. 

Both SOC approximation and the fact that elementary diagrams are neglected may lead to a 
violation of the variational principle. To keep this effect under control, an additional constraint 
on the kinetic energy with respect to the optimization performed in |24] has been considered. 




Figure 4. Energy per particle for PNM, obtained using the density-dependent potential added 
to the Argonne v'^ (a) and to Argonne Vq (b) potentials. The energies are compared to those 
obtained from the genuine three-body potential and from the two-body potentials alone. 

We carried out AFDMC simulations [T] for PNM with A = 66 and SNM with ^ = 28 nucleons 
in a periodic box. The finite-size errors in PNM simulations have been investigated in [25] by 
comparing the Twist Averaged Boundary Conditions (TABC) with the Periodic Box Condition 
(PBC). The kinetic energy of 66 fermions approaches the thermodynamic limit very well, so 
that the energies of 66 neutrons computed using either TABC or PBC turn out to be almost 
the same. We can infer that the finite-size errors in the present AFDMC calculations for PNM 
do not exceed 2% of the asymptotic value of the energy. 

The finite-size effects of SNM calculations can be estimated from the difference of the energies 
of PNM obtained with 14 neutrons and the TABC asymptotic value, which is of the order of 
7%. 

FHNC/SOC and Monte Carlo calculations provide very close results, as shown in Figs. H] 
and [5l to be compared with those of Ref. [2] where the agreement between FHNC and Monte 
Carlo methods were not nearly as good. 

The EoS of PNM, displayed in Fig. SI obtained with the three-body potential UIX and 
using the density-dependent two-body potential are very close to each other. In the same figure, 



for the sake of comparison, we also report the results of calculations carried out including the 
two-body potential only. 

Despite cluster contributions proportional to have been neglected in the density dependent 
potential, with the exception of the line with diamonds of Fig. [21 even at high densities the EoS 
relative to vi2{p) remain close to those obtained with UIX. Probably, in this case a compensation 
among second and higher order terms takes place. 




Figure 5. Same as in Fig. [H but for SNM 



The density-dependent potential has been also employed in AFDMC calculations. As can be 
seen in Fig. U the triangles representing the results of this calculation are very close, when not 
superimposed, to the circles corresponding to the UIX three-body potential AFDMC results. 

For what concern the EoS of symmetric nuclear matter, see Fig. [5l at densities lower than 
p = 0.32 fm~^, the curves resulting from UIX and the density-dependent potential are very close 
to one other, while for p > 0.32 fm~^ a gap between them appears. 



Table 1. Values for the saturation densities, the binding energy per particle, and the 
compressibility of SNM relative to the EoS of Fig. [5j 



FHNC/SOC 


v'q + Ui23 


v'q + v{p) 


+ ^123 


4 + v{p) 


AFDMC 


+ v{p) 


Po (fm-^) 


0.17 


0.16 


0.16 


0.15 


Po (fm-^) 


0.17 


Eq (MeV) 


-11.3 


-11.2 


-10.3 


-10.3 


Eo (MeV) 


-10.9 


K (MeV) 


205 


192 


189 


198 


K (MeV) 


201 



In Table [H the saturation density pQ, the binding energy per particle E{pq) and the 
compressibility K = 9po{dE(p)/dp)'^ for all the EoS of Fig. [5] are listed. It is remarkable that 
the values obtained with density dependent potential are very close to those resulting from the 
genuine UIX three body potential. While the saturation density is well reproduced, this being 
not surprising since the constant Uq of the UIX potential has been fitted to this value, the values 
of the compressibility are slightly lower than the experimental ones. AFDMC calculations lower 
the variational results for Eq by 0.6 MeV only, showing that the binding energy overstimation 
of the FHNC/SOC calculations performed with UIX potential can not be solved with a better 
choice of the variational wavefunction. 



6. Conclusions 

We have developed a novel scheme, suitable to obtain an effective density-dependent NN 
potential taking into account the effects of three- nucleon interactions. Our approach is fully 
consistent with the treatment of correlations underlying the FHNC and AFDMC approaches. 

The PNM and SNM equation of state resulting from the density-dependent potential turn out 
to be very close to those obtained with the UIX three-body potential. In this context, a critical 
role is played by the treatment of both dynamical and statistical correlations, a distinctive feature 
of our approach, as compared to different reduction schemes based on effective interactions 

For the first time, a AFDMC calculation of the equation of state of SNM consistently including 
the effects of three nucleon forces has been carried out. The results of this calculation show 
that the fg-|-UIX hamiltonian, or equivalently the one including the effective potential, fails 
to reproduce the empirical data. The discrepancy has most likely to be ascribed either to 
deficiencies of the UIX model or to the effect of interactions involving more than three nucleons. 

As a further development, we are studying the dependence on the specific model of three- 
nucleon force [30|, as well as the inclusion of four- and many- nucleon interactions, whose effects 
are expected to be critical for the determination of the properties of high density neutron star 
matter . 

Our density dependent potential could be easily employed in many-body approaches other 
than those based on the CBF formalism or quantum Monte Carlo simulations, such as the 
G-matrix and self-consistent Green function theories |28[ [29l [3T] . 



References 

[I] K E Schmidt and S Fantoni, Phys. Lett. B446, 99 (1999). 

[2] S Gandolfi, F Pederiva, S Fantoni and K E Schmidt Phys. Rev. Lett. 98, 102503 (2007). 

[3] A Sarsa, S Fantoni, K E Schmidt F and Pederiva , Phys. Rev. C 68, 024308 (2003). 

[4] I Bombaci, A Fabrocini, A Polls and I Vidana, Phys. Lett. B 609, 232 (2005). 

[5] B S Pudliner, V R Pandharipande, J Carlson and R B Wiringa, Phys. Rev. Lett. 74, 4396 (1995). 

[6] S C Pieper, Rivista del Nuovo Cimento 31, 709 (2008). 

[7] I Lagaris and V R Pandharipande, Nucl. Phys. A359, 349 (1981). 

[8] B Friedman and V R Pandharipande, Nucl. Phys. A361, 502 (1981). 

[9] A Lovato, O Benhar, S Fantoni, A Yu Dlarionov, and K E Schmidt, Phys. Rev. C 83, 054003 (2011). 
[10] A Kievsky, M Viviani, L Girlanda, and L E Marcucci, Phys. Rev. C 81 044003 (2010). 

[II] O Benhar, V R Pandharipande and S C Pieper, Rev. Mod. Phys. 65, 817 (1993). 

[12] V R Pandharipande, I Sick, P K A. DeWitt Huberts, Rev. Mod. Phys. 69, 961 (1997). 

[13] S Fantoni and S Rosati, Nuovo Cimento A20, 179 (1974). 

[14] O Benhar and M VaUi, Phys. Rev. Lett. 99, 232501 (2007). 

[15] O Benhar, N Farina, S Fiorilla and M Valh, AIP Conf. Proc. 1056, 248 (2008). 

[16] B S Pudliner, V R Pandharipande, J Carlson, S C Pieper and R B Wiringa, Phys. Rev. C 56, 1720 (1997). 

[17] R B Wiringa, V G J Stokes, R Schiavilla, Phys. Rev. C 51, 38 (1995). 

[18] R B Wiringa, S C Pieper, Phys. Rev. Lett. 89, 182501 (2002) 

[19] J Fujita and H Miyazawa, Prog. Theor. Phys. 17, 360 (1957). 

[20] J W Clark, Prog. Part. Nucl. Phys. 2, 89 (1979). 

[21] S Fantoni and V R Pandharipande, Phys. Rev. C37, 37 (1988). 

[22] V R Pandharipande and R B Wiringa, Rev. Mod. Phys. 51, 821 (1981). 

[23] S Kirkpatrick, C D Gelatt and M P Vecchi, Science New Series, 220, 671 (1983). 

[24] R B Wiringa, V Fiks, and A Fabrocini, Phys. Rev. C 38, 1010 (1988). 

[25] S Gandolfi, A Yu Illarionov, K E Schmidt, F Pederiva and S Fantoni, Phys. Rev. C79, 054005 (2009). 

[26] K Hebelerl and A Schwenk, Phys. Rev. C 82, 014314 (2010). 

[27] J W Holt, N Kaiser and W Weise, Phys. Rev. C 81, 024002 (2010). 

[28] M Baldo e F Burgio, in Physics of Neutron Stars Interior. Eds. D Blaschke, N. K. Glendenning and A. 

Sedrakian. Lecture Notes in Physics 578, 1 (2001). 
[29] A. Rios, A. Polls and I. Vidaiia, Phys. Rev. C 79, 025802 (2009). 
[30] A. Lovato, O. Benhar, S. Fantoni, and K. E. Schmidt, in preparation. 
[31] W.H. Dickhoffand C. Barbieri, Prog. Part. Nucl. Phys. 52, 377 (2004). 



